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We study massless and massive graviton modes that bind on thick branes which are supergravity 
domain walls solutions in D- dimensional supergravity theories where only the supergravity multiplet 
and the scalar supermultiplet are turned on. The domain walls are bulk solutions provided by 
tachyon potentials. Such domain walls are regarded as BPS branes of one lower dimension that are 
formed due to tachyon potentials on a non-BPS D-brane. 



m 
O 
O 
(N 

< 



> 

oo 
O 

m 
o 



Oh 



X 



PACS numbers: ll.10.Kk, 04.50.+h, ll.25.-w, 11.27.+d 



I. INTRODUCTION 

The localization of gravity on brane worlds is a very 
active matter of investigation in modern physics that 
was initiated in 1]. We study here background of p- 
branes created by scalar fields of the bosonic sector of 
a D-dimcnsional supergravity theory. We shall mainly 
explore the localization of the graviton on thick p-branes 
0, 01 supported by these scalars. These thick p-branes 
are domain wall solutions extended in p = D — 2 dimen- 
sions. 

As one knows, thick branes can be studied with- 
out adding any scalar field, by smoothing out singular 
brane solutions However, we are interested in a su- 
pergravity theory with the contents relevant to realize 
brane/string cosmology in arbitrary D-dimensions where 
scalar fields can play important role. We recall that such 
scalars play the role of tachyon fields whose dynamics are 
truncated up to first order on the derivatives. We choose 
some particular potentials in order to have connection 
with modern cosmology. One of the many issues found 
in cosmology is the accelerated expansion of our universe 
0, and this may have to do with dark energy. In view 
of this problem new proposals of cosmological scenarios 
as, for instance, quintessence jfj and tachyon cosmology 
0,0 HI have been considered in the literature. 

In the present paper we focus mainly on the scalar con- 
tents of tachyon cosmology in a supergravity theory. We 
consider here two different tachyon potentials which we 
call type I and type II. The tachyon potential of type I 
that one uses here is a potential in the bulk, for which 
the global minima are at infinity [Tol . The tachyon 
potential of type II is unbounded from below but there 
are BPS solutions, i.e., domain wall solutions, with fi- 
nite energy, connecting inflection points of the potential 
173 . These inflection points defines an interval in which 
the potential engenders the properties of being positive 
definite, with global minima at the boundaries of the 
interval. The tachyon potentials of type I and type II 
that we study are connected by a deformation function 
[l2j . They correspond to different and relevant cosmolog- 
ical scenarios. For models which minima are at infinity 
the tachyon matter is pressureless and may be consid- 



ered a cold dark matter candidate 01 On the other- 
hand, for potentials that present minima at finite criti- 
cal points the tachyon matter has negative pressure and 
may be considered a candidate for quintessence 0. Fi- 
nally we should mention that in spite of the fact that the 
tachyon potential of type II considered here in general 
breaks the supersymmetry, one sees that one half of the 
supersymmetry is preserved in the solitonic sector, i.e., 
in the sector where fields configurations are in the finite 
interval and also where the domain walls (thick p-branes) 
can form. 

The paper is organized as follows. In Sec. ^we intro- 
duce the D-dimensional supergravity theory. In Sec. IIIII 
we discuss the issue of graviton mode localization on p- 
branes. The properties of the graviton localization for 
the backgrounds created by the aforementioned poten- 
tials are discussed in Sec. HVI and in Sec. El We end the 
paper in Sec. IVII where we present our final considera- 
tions. 



II. SUPERGRAVITY ACTION IN ARBITRARY 
D-DIMENSIONS 

We consider a supergravity Lagrangian for arbitrary 
spacetimes in _D-dimensions (D > 3) invariant under su- 
persymmetry transformations, at least to lowest order in 
the fermions, of the form [lj, []J| 
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where k = 1/M* is the D-dimensional Planck length, 
A/* is the fundamental D-dimensional Planck scale, e = 
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dete^j = | det gMN] 1 ^ 2 , with a mostly plus signature 
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W is the superpotential, the pair of graviton and grav- 
itino fields (e^ jV'm) corresponds to the supergravity mul- 
tiplet whereas the pair of scalar and spin- 1/2 fermion 
fields {4>, x l ) forms a scalar supermultiplet. We use Tm to 
represent Dirac matrices. We include an internal spinor 
index i for generality. The action Q is invariant under 
the following supersymmetry transformations: 
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where e l is a local supersymmetry parameter and 
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A similar consideration but exploring other issues con- 
cerning localization of gravity on brane in supergravity 
has also been considered in [T^, [T3, 13 • 

For our purposes we employ the "generalized" Randall- 
Sundrum metric: 



ds 
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e Aiv) r)„ v dx^dx v + dy 2 , 



(10) 



where fi,u = 0, 1, 2, D — 2 are indices on the (D — 2)- 
brane. For this metric the Killing spinor equations 
8ij) % M = and 8\ l = provide us the following equa- 
tions: 



d y A = t2k D ~ 2 W, 
d y <j> = ±2(D-2) 



dW 



(11) 
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Notice that we are assuming the scalar field and the 
spinor field e l only depend on the transverse coordinate y. 
These equations describe BPS solutions since they come 
from the fermion supersymmetry transformations. Such 
solutions preserve only half of the supersymmetries. 



III. GRAVITON MODES ON (D 2)-BRANES 



where $ describe the wave function of the graviton on 
non-compact coordinates. M, N — 0,1, 2, D — 1. This 
can be found by taking into account only the trace- 
less transverse (TT) sector of the linear perturbation, 
in which the gravity equations of motion do not couple 
to matter fields. 

Let us consider $ = h(y)M(x fJ ') into (|13fl and the fact 
that D D _iM — m 2 M, where D D _iM is the flat Lapla- 
cian on the tangent frame. Thus the wave equation for 
the graviton through the transverse coordinate y reads 



d y (V=99 vy dyHy)) 



-m 2 \g 00 \h(y). 
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This is our starting point to investigate both zero and 
massive graviton modes on the branes. 

Using the components of the metric <|10(l into the equa- 
tion (|14|) we have 

l -{D-l)d v Ad y h(y) + d 2 h(y) = -m 2 e- A ^h(y), (15) 

which can be written as a Schrodinger like equation by 
changing the metric (|1UH to a conformally flat metric as 



ds 2 D = e^Xri^dx^dx" + dz 2 ), 



(16) 



and also by considering the following changes of variables: 
h(y) = V(^)e" A(z)( - D - 2)/4 and z(y) = J e~ A ^/ 2 dy. In this 
way we can write l|15f) as a Schrodinger like equation 

- d 2 ip{z) + V(z)ip(z) = m 2 iP(z) (17) 

where the potential V(z) is given in the general form 

D-2 



V(z) = (d z A) 2 
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Here we notice that the potential depends only on the 
geometrical variable A(z) up to constant factors. It is 
easy to check that for D = 5 and A(y) — —2k\y\ one 
reproduces the Randall-Sundrum Q] potential. 



IV. EXPLICIT EXAMPLES OF TACHYON 
POTENTIALS 

We shall look for soliton solutions provided by poten- 
tials of tachyons that live on unstable D-branes in su- 
perstring and bosonic string theory 0, I22I |23| — for 
D-brane/anti-D-brane systems, the tachyon field should 
be complex, and this will not be considered here. The 
tachyon dynamics on the world-volume of a p-brane is 
described by the action 0, |25| 



The wave function of graviton modes due to linearized 
gravity equation of motion in an arbitrary number of 
dimensions (D > 3) is given by 0, 0, |2l| 
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where V(T) is the tachyon potential, which should have 
the general property of being a positive definite function 
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with stable minima at \T\ — > ±00. Expanding this action 
up to quadratic first derivative terms gives [26j 

dP +1 xV{T)(\ + ±Tf v d ll Td v T+ ...) 

-\^ v d ll 4>d v 4>-v{cj>) + ..), (20) 

where in the second equation above we assumed that 
V(T) = 4>t — nere 't'T stands for dr4> — such that we 
can replace fad^T to d^4>. The tachyon dynamics is now 
described by the action for the scalar field <fi. This action 
can be viewed as the scalar sector of the action {Q . The 
scalar field here is related to the former tachyon field T 
and transforms as (f> = f(T). 

We now focus on two special examples that allow to 
form domain walls solutions. 

The model type I: <fi — (1/A) arcsinh (AT). This pro- 
duces a tachyon potential 



1 



1 + A 2 T 2 



= sech 2 (A</>). 
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Notice that the potential above has the property of hav- 
ing global minima at T(or <fi) = ±00. Potentials of this 
type are related to the massless limit of supersymmet- 
ric QCD. In this limit, supersymmetric QCD coupled to 
supergravity might play the role of the hidden sector of 
supergravity theories [23]. These potentials have been 



employed in models of cosmology |28 

The model type II: <\> = T/y/1 + T 2 . This produces a 
tachyon potential of the form 



V{T{<j>)) = 



1 
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We see here a different asymptotic behavior, since the 
above potential is unbounded below. However, the two 
inflection points at = ±1 identify the interval — 1 < <f> < 
1, where BPS solution exists [lj- Thus, if one restricts 
the field to the interval 6 [—1,1], the potential gets the 
required profile. 

The above models are relevant to cosmology, where 
each one plays a different role 0. For the purposes of 
this paper, we are mainly concerned with the fact that 
these potentials have "minima" at finite or infinite crit- 
ical points. They allow for domain walls solutions that 
can be reg arded as BPS branes of one lower dimension 
— see [23] and references therein. 

Now we shall look for domain walls solutions and their 
gravitational localization properties. As a first example 
we shall consider the model type I. We introduce a super- 
potential with which we are able to find soliton solutions 
that connect the superpotential critical points at infinity. 
Let us consider the general superpotential 
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7j arctan[sinh(A0)]. 
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By using © we can see that (|23|l produces the potential 
1/(0) = 40D-2) 2 {^J sech 2 (A0) 
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■ arctan 2 [sinh(A</>)] j . (24) 
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The first term in l|24|l comes from the derivative 
— = -sech(A0) 



(25) 



and should appear even in global supersymmetry jllj . 
i.e., when k = 0. The potential (|24|) (for k = 0) up to 
constant factors is of the type given in 121|) . 

We can deal with the system of equations (jTT)l and (fT2|) 
easily because they are decoupled. The soliton solution 
can be easily found by using ltT2l : 



<t>(0 = ±-arcsinh£, 



(26) 



where we have considered £ = 2(D — 2)//j/. These solu- 
tions asymptote to infinity. Now we apply the positive 
solution to the equation that describes the geometry (jl 1|> 
— we have set here n = 1 for simplicity. By using the su- 
perpotential (|23|l we find 



A(0 =2[m(l + £ 



2x1/2 



£ arctan £ 



(27) 



We set 2A 2 (Z3 — 2) = 1 It is not difficulty to see that the 
warp factor e A ^ falls off to zero as |£| goes to infinity. 

We now search for the conformal "z-coordinate" in 
order for to write down the potential (|18fl in the 
Schrodinger like equation (fTTjl . We need to perform the 
integral z(y) — J e~ A ( v ^ 2 dy, i.e., 



z{y) = 



2{D 



[{ arctan£-ln (1+£ 2 ) 1/2 ] 



(28) 



For very large [A, the thin wall limit, the integral in 
(|28|l turns out to have approximately the simple form 
/ exp = sgn(£) exp (]£]). By continuity in the pa- 

rameter fi, this means in the thick wall regime we can 
approximate such integration by the function sinh£, by 
smoothing out back the functions sgn(£) — > tanh£ and 
|£| — >ln cosh£. 

All we have assumed above suffices to guarantee in- 
tegr ability and inversion in J33J). In this sense, up to a 
constant we have set to zero, we have the following: 

, , sinh lauv] , . arcsinh \auz] 
g ( y ) = Li^i or y ( z ) = U_J, (29) 



07/ 



where a = 2(D — 2). The approximation (|29|l can actually 
be confirmed numerically, since the numerical integral of 
has a "sinh" type profile. 
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The potential (|18(l is now written explicitly as 
1 a^/j 2 arctan 2 [arcsinh(a^z)] 



V(z) 



16 1 + a 2 fj?z 2 

1 a /j, 3 arctan[arcsinh(a/xz)]z 
4~ 
1 



(1 + aVz 2 ) 3 / 2 



4 (1 + a 2 /^ 2 z 2 )[arcsinh (aifiz) + 1] 



(30) 



This potential has the asymptotic behavior: V(z = 
±oo) = 0. This means that the potential provides no 
mass gap to separate the graviton zero mode from KK 
(massive) modes. At z = the potential has a minimum 
with depth V(0) = — (l/4)a 3 /i 2 . In fact, this is a volcano 
type potential, which tends to the singular one found in 
the Randall-Sundrum scenario as /j,—>oo 

The analysis for the tachyon potential of the model 
type II is completely analogous to the previous case. The 
superpotential is chosen to be of the form 
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which is restricted to the interval —a < <f> < a and pro- 
duces the potential 
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If gravitational effects are not included in the theory, i.e. 
k = 0, the remaining term in (|32H . which has also been 
found in global supersymmetry 12], is of the type in 
the Ea. (|22|l . This term is in fact related to the remain- 
ing term of the potential 1(24(1 . for n = 0, by applying a 
function of deformation f((f>) = a tanh (A</>) as has been 
shown in |l2j . Precisely, we mean that given the tachyon 
potential of the model type II V K=0 (4>) ~ la A )(a 2 — 2 ) 3 
one can obtain the tachyon potential of the model type 
I V k= q(4>) ~ (f^/X 2 ) sech 2 (A0) by just using the transfor- 
mation 



V K=0 {4>) = 



v K = [f(4>)} 
\f'{4>)] 2 



(33) 



As we mentioned in the introduction, the tachyon po- 
tential of the model type II preserves half of supersymme- 
tries in the solitonic sector, where there are BPS solutions 
defined in the interval —a<4> <a. 

We use the equation i|12[l and the tachyon superpo- 
tential of the model type II (|31|l to obtain the following 



solution 



= ±o- 



i 



(34) 



Recall that we are considering £ = 2(D — 2)fj,y. This is a 
topological soliton solution connecting two critical points 
of the superpotential W which are inflections points of 
V K= o(0). This solution, just as the domain wall solution 
found in the previous case, is regarded here as a "thick" p- 
brane, where y represents the only transverse coordinate 
in D=p+ 2 dimensions. 

We use equation (|ll(l , setting k = 1 for simplicity, in 
order to find the warp factor on this background. We get 



A(0 - 2 



1/3 



£ arctan £ 




,(35) 



where we have set a = 4^/(1? - 2)/3. Notice that e A(?) 
falls off to zero as |£| goes to infinity. The same arguments 
used before can be applied to find that here the conformal 
z-coordinate is also given by the equation (|2T)|) . 

The potential (jlSJ) can now be written explicitly as 



[ ' 4 F 2 



16 



2F 4 F$z 



3 (F 2 + 1 



?F 2 



, (36) 



where F\ =arcsinh[2(_D — 2)(iz], F3 = D — 2 and 
[l + 4(£-2)Vz 2 ] 1 /2 
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F 2 

Fi = -2 arctan (Fi) - 



2Fi(3F 1 2 + l) 
3 (F 2 + l) 2 ' 



This potential as the potential 1)30(1 has the asymptotic 
behavior: V(z — ±00) = 0. This means that there is 
no mass gap. At z — the potential has a minimum 
with depth V(0) = -(16/3)(D - 2) 3 /i 2 = -(2/3)a 3 ^ 2 ~ 
compare this with the depth in the previous case. In the 
following we shall look for solutions which have either 
zero or massive eigenvalues. 

The graviton zero mode on the (D — 2)-brane can be 
found easily since the Schrodinger equation 1(17(1 with 
the general potential ((18(1 can be written as Hip = m 2 t/j, 
where the Hamiltonian operator is given by H = Q]Q, 
with Q = -d z + (1/4) (D - 2)d z A — see || H |H |Hl 
for further details on this issue. Since the operator H is 
positive definite, there are no normalizable modes with 
negative energy (mass), i.e., there is no tachyonic gravi- 
ton mode. Thus the stability of the (D — 2)-brane is en- 
sured and the graviton zero mode is the lowest mode in 
the spectrum. The operator Q annihilates the zero mode 
ipo(z), i.e. Qipo{z) = 0. This implies that the graviton 
zero mode is 



V>o(z) 



N a e— A{ - Z \ 



(37) 



where Nq is a normalization constant. Notice that the 
metric warp factor e A ^ for both cases studied above 
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asymptotes to zero, even in the conformal z-coordinate, 
as \z \ — >oo, such that the graviton zero mode l|37[) is nor- 
malizable, i.e. \^o{z)\ 2 dz < oo. We then conclude 

that the graviton zero mode binds to the (D — 2)-brane 
located at z — 0. In addition to such a mode, the po- 
tentials l|3U|) and H36|) suggest the existence of massive 
modes (KK modes). These modes can affect the New- 
tonian gravitational potential between two masses at the 
(D— 2)-brane located at z — 0. In the next section we shall 
look for correction to the Newtonian potential. Due to 
the difficulty of integrability in (|17fl for massive modes, 
the calculations that follow involve only an asymptotic 
analysis. 



V. CORRECTION TO THE NEWTONIAN 
POTENTIAL AT THE 3-BRANE 

The four-dimensional gravitational potential is cor- 
rected by massive Kaluza-Klein modes 1]. Since the KK 
modes correct the Newtonian potential, the correction 
should be very small, at least for the case of phenomeno- 
logical interest, i.e. D — 5, in order to keep effectively the 
well known Newtonian potential in our four-dimensional 
universe, the 3-brane. 

We can estimate the static Newtonian potential be- 
tween two particles of mass mi and m,2 on the 3-brane 
by considering the exchange of the zero-mode and con- 
tinuum KK mode propagators Q. The general form of 
this potential in five dimensions |2J, |j] is given by 



V(r) ~G 



77117712 



-M. 



N ' 



_ 3 77117772 



7" 



dm t 



> m (0)| 2 , (38) 



where M* is the fundamental Planck scale in D — 5 di- 
mensions. Notice we are integrating on the masses of the 
continuum modes. The lower limit mo in the integral is 
set to zero when there is no mass gap separating such 
continuum modes from the zero mode. 

In Ea. H38|l we see that in order to compute the cor- 
rection of the Newtonian potential we need to specify 
the ip m (z) functions for all possible m. In the models we 
are studying these functions cannot be found analytically 
in a closed form, therefore we only consider an asymp- 
totic analysis of the problem. In this sense, we study 
the quantum mechanics problem for the potentials l|3U|) 
and H3tjf) . by looking for solutions in the asymptotic limit 
\z \ 3> 1/afi. In such a limit the potentials l|3U|) and i|3t)|) 
asymptote to 



V(z) 



a'(a' + 1) 



(39) 



where a' — an/8. We can now solve the quantum me- 
chanics problem for these potentials |4j. The solution is 
given by a linear combination of Bessel functions: 

ip m {z) = a m z 1/2 Y a , +1 / 2 (mz) + b m z 1/2 J a , +1/2 {mz).(A0) 



After considering appropriate boundary conditions it is 
possible to estimate the coefficients a m and b m , such that 
one can estimate the value of the wavefunction at z = 0. 
The leading term in m is given by 



*.«•>- (t) 



a'-l 



(41) 



The factor of k should agree with dimensional analysis. 
We can here approximate this factor to the inverse of 
width of the domain walls, i.e., which has dimen- 

sion of mass. 

We can now apply this result to Eq.(j3Hl to find the 
correction to the Newtonian potential. Integrating over 
all masses by using (|41() . we find the four-dimensional 
gravitational potential corrected as 



V(r) = G 



?71l?7l2 



JV " 



7" 



C 



(At) 



2d/ 



(42) 



where C is a dimensionless number and Gjv ~ k/M^. 
Notice that the correction is highly suppressed for k 
around the fundamental five dimensional Planck scale 
and r around the size tested with gravity []]. Since 
a' = air/8, recalling that a — 2(D — 2) we have the follow- 
ing situations: For the model type I, recall we have set 
2X 2 (D - 2) = 1, thus a' = 7r/8A 2 . On the other hand, for 
the model type II, recall we have set a = 4^/ (d — 2)/3, 
thus a' — (3/64)7ra 2 . At large distances, the Newtonian 
potential is not modified for a' > 1/2. This implies the 
limits < A < \px /2 and a > (32/37T) 1 / 2 on the parame- 
ters. 

Notice that in both cases discussed above, the inverse 
power law of the correction is not always like l/r™, where 
n is an integer number. For arbitrary choice of both pa- 
rameters A and a we can have rational inverse power law. 
In this situation we are led to think of "fractal dimen- 
sions" inducing correction to the Newtonian potential. 
For instance if A = (l/6)v / 37r, a' — 3/2, as in the Randall- 
Sundrum case , then the correction goes as l/r 2 . How- 
ever, for A slightly larger than this value, we have l/r p , 
where 1 < p < 2. On the other hand, for a — 4^2/tt we 
get the power law 1 /r 2 , whereas for a slightly larger than 
this value, we find l/r p , with 2<tj<3. 



VI. CONCLUSIONS 

In this paper we have investigated the localization of 
zero and massive graviton modes on thick branes sup- 
ported by tachyon fields that live on the world-volume 
of non-BPS D-branes. The decay process of such ob- 
jects due to tachyons ends up creating BPS branes of one 
lower dimension via tachyon fields — see (2^ and refer- 
ences therein. If we consider that the process happens 
in D = 5 dimensions, then the tachyon dynamics in the 
action can be regarded as the tachyon dynamics on 
the world-volume of a non-BPS D4-brane of superstring 
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theory. The brane of one lower dimension here is a BPS 
3-brane and might be our four-dimensional universe. 

In the quantum mechanics problem for the graviton 
spectrum on the BPS 3-brane, we have found volcano 
potentials for both models type I and type II. These po- 
tentials have just one bound state, the zero mode, and a 
continuum of Kaluza-Klein modes that correct the New- 
tonian potential in a highly suppressed way. With re- 
spect to the matter field excitations on the p-branes, 
it has already been shown llj that there are no mat- 
ter zero modes on p-branes solutions in the supergravity 
theory that we have considered. This is because there 
are pre-factors like e tyD ~ n ^ A ^W{z), for n = 1,3, in the 
(D — l)-dimensional effective action that come from the 
integration of the matter zero modes. Such factors are 



zero due to exponential fall-off of the warp factor e A< - z \ 
The investigation of matter massive modes is out of the 
scope of the paper. Such modes might reveal some new 
interesting physics and we leave this to be explored else- 
where. 
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